Constraining the growth factor with baryon oscillations 
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The growth factor of linear fluctuations is probably one of the least known quantity in observa- 
tional cosmology. Here we discuss the contraints that baryon oscillations in galaxy power spectra 
from future surveys can put on a conveniently parametrized growth factor. We find that spectro- 
scopic surveys of 5000 deg^ extending to z « 3 could estimate the growth index 7 within 0.06; a 
similar photometric survey would give A7 ~ 0.15. This test provides an important consistency 
check for the standard cosmological model and could constrain modified gravity models. We discuss 
the errors and the figure of merit for various combinations of redshift errors and survey sizes. 



I. INTRODUCTION 

The characterization of dark energy (DE) has been so 
far based almost uniquely on background tests at rather 
low redshifts {z < 1.5: Riess et al. 1998, Perlmutter et 
al. 1999, Tonry et al. 2003, Riess et al. 2004, Astier et 
al. 2006, Eisenstein et al. 2005) or very large redshifts 
{z « 1000: e.g. Netterfield et al. 2002, Halverson et al. 
2002, Lee et al. 2002, Bennett et al. 2003, Spergel et al. 
2006). These tests are based essentially on estimations 
of luminosity Dl (z) or angular-diameter distances D{z), 
i.e. on integrals of the Hubble function H{z) which, in 
turn, contain integrals of the equation of state. Only very 
recently tests involving the linear perturbations have be- 
gun to be discussed, using methods based on the inte- 
grated Sachs- Wolfe effect, weak lensing and high-redshift 
power spectra (e.g. Boughn & Crittenden 2004, Refregier 
et al. 2006, Crotts et al. 2005 ). However, it is fair to say 
that the growth function is still one of the least known 
quantity in cosmology. So far, it is possible to quote only 
two pubHshed results that put Hmits on it: the value at 
z « 0.15 obtained in 2dF (Hawkins et al. 2003; Verde 
et al. 2002) and the z « 3 result from Lyman— a clouds 
(McDonald et al. 2005). Defining G{z) = S{z)/S{0) {S 
being the matter density contrast) we have for 



/ 



dlogG 
dloga 



(1) 



the value / = 0.51 ± 0.15 for 2dF at z « 0.15 and 
/ = 1.46 ± 0.29 for the Lyman-a at z « 3 . These re- 
sults show clearly how large is the degree of uncertainty. 
Actually the uncertainty is much larger than it appears 
from the quoted statistical errors. In the case of the low- 
z estimate, the result is obtained by estimating the bias 
from higher-order statistics, which is known to be partic- 
ularly sensitive to the selection effects, to incompleteness 
etc.; different methods give in fact quite different results 
(see discussion in Hawkins et al. 2003). In the case of 
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the high-z estimation, the main problem is the recon- 
struction of the bias factor from numerical simulations 
which, by their nature, are performed only in a limited 
range of fiducial models. It is therefore important to test 
the growth factor with other methods and with improved 
datasets. 



A test of the growth factor would be important both as 
a consistency check for the standard cosmological model 
(since / is determined by H{z) in a standard cosmol- 
ogy) and as a constraint on non-standard models like 
e.g. modified gravity. In fact, models that modify the 
Poisson equation will also generically modify the pertur- 
bation equation for the matter density contrast S. As an 
example, models in which dark energy is coupled to mat- 
ter display a growth index which deviates from the stan- 
dard case at all epochs (see e.g. Amendola & Tocchini- 
Valentini 2003; Demianski et al. 2004; Nunes & Mota 
2004). Several other papers discussed the parametriza- 
tion of the perturbation equations in modified gravity 
models, see e g. Ishak et al. (2005), Heavens, Kitch- 
ing and Taylor (2006), Taylor et al. (2007); Heavens, 
Kitching, Verde (2007) , Caldwell, Cooray and Melchiorri 
(2007), Amendola, Kunz and Sapone (2007), Zhang et al. 
(2007). 



In this paper we investigate the extent to which baryon 
oscillations can set limits to G(z) in future large-scale ob- 
servations at z up to 3. The method we use is based on 
recent proposals (Linder 2003, Blake & Glazebrook 2003, 
Seo & Eisenstein 2003) to exploit the baryon acoustic os- 
cillations (BAOs) in the power spectrum as a standard 
ruler caHbrated through CMB acoustic peaks. In partic- 
ular, Seo & Eisenstein (2003; SE) have shown the feasi- 
bility of large (100 to 1000 square degrees) spectroscopic 
surveys at z « 1 and z w 3 to put stringent limits to 
the equation of state w{z) and its derivative. As it is 
well-known, BAOs have been detected at low z in SDSS 
(Eisenstein et al. 2005); the detection at large z , where 
more peaks at smaller scales can be obtained, is likely to 
become one of the most interesting astrophysical endeav- 
ours of the next years. 
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II. BACKGROUND EQUATION 

Here we review the basic equations and notation for 
the background evolution and for the Unear fluctuations. 
The evolution of the dark energy can be expressed by the 
present dark energy density CIde and by a time- varying 
equation of state (see Copeland, Sami & Tsujikawa 2006 
for a recent review): 



wiz) = - 
P 



(2) 



Given w{z), the dark energy density equation is p{z) = 
p(0)a-3(i+*) where 



w{z) 



log(l 



i + z)Jol + z' 



The Hubble parameter and the angular diameter dis- 
tance, H{z) and Da{z), assuming a flat universe fi^ + 
CIde = 1, become respectively: 

(z) = H^[n^,{l + zf + (1 - + (4) 

and 

dz' 



Da{z) 



H{z') 



where the total matter density is 



3u) 



(5) 



(6) 



It is well known that a good approximation to the 
growth index for sub-horizon scales in flat models is given 
by (Lahav et al. 1991, Wang and Steinhardt 1998 ) 



_ aiogG ^ 

J = ^1 = {a) 

a log a 



(7) 



This introduces a new parameter 7, beside those that 
characterize the background model (see also Linder 2005, 
Percival 2005) . 

We remark that a recent analysis of most of the ex- 
tant data produced the result 7 = 0.6lo!3 (Di Porto & 
Amendola 2007) . 



III. FISHER MATRIX FORMALISM 

Following Seo & Eisenstein (2003; hereinafter SE) we 
write schematically the observed galaxy power spectrum 
as: 



D\{z)Hr{z) 

+Pshot{z) 



G^{z)b{zy{l + PnYPor{k) 

(8) 



where the subscript r refers to the values assumed for the 
reference cosmological model, i.e. the model at which we 



evaluate the Fisher matrix. Here Pshot is the shot noise 
due to discreteness in the survey, /z is the direction cosine 
within the survey, Pq is the present spectrum for the 
flducial cosmology. For the linear matter power spectrum 
we adopt the flt by Eisenstein & Hu (1999) (with no 
massive neutrinos and also neglecting any change of the 
shape of the spectrum from small deviation around w = 

-!)• 

The wavenumbcr k is also to be transformed between 
the fiducial cosmology and the general one (SE; see also 
Amendola, Quercellini, Giallongo 2004, hereinafter AQG, 
for more details). The bias factor is defined as: 



b{z) 



^711 (^) 



(9) 



and for the fiducial model is estimated by comparing the 
8Mpc//i cell variance as,g of the galaxies corrected for 
the linear redshift distortion with the same quantity for 
the total matter. Clearly, the growth function is degener- 
ate with the bias except for the redshift correction factor 
{l+Pn"^). Since we marginalize over (3, it is clear that the 
redshift correction plays a crucial role for as concern the 
estimation of the growth factor. The linear correction we 
use should therefore be considered only a first approxi- 
mation and more work to go beyond Kaiser's small-angle 
and Gaussian approximation is needed, as discussed in 
Hamilton & Culhanc (1996), Zaroubi & Hoffman (1996), 
Tegmark et al. (2004) and Scoccimarro (2004). 

The total galaxy power spectrum including the errors 
on redshift can be written as (SE) 



P {z, k) = Pobs {z, k) e 



,222 



(10) 



where cxr 



Sz 
HIT) 



is the absolute error on the measure- 



ment of the distance and 6z is the absolute error on red- 
shift. Given the uncertainties of our observations, we 
now want to propagate these errors to compute the con- 
straints on cosmological parameters. The Fisher matrix 
provides a useful method for doing this. Assuming the 
likelihood function to be Gaussian, the Fisher matrix is 
(Eisenstein, Hu & Tegmark 1998; Tegmark 1997) 



2tt 



rkmaa 



dl0gP{kn)dl0gP{kn) 



■V, 



eff ■ 



87r3 



■dk 
(11) 

where the derivatives are evaluated at the parameter val- 
ues of the fiducial model and 14/ / is the effective volume 
of the survey, given by: 



eff 



n{f)P{k,n) + 1 
n (r) P (fc, 12) 
n(r)P(fc,M) + l. 



df = 



survey 



(12) 



where the last equality holds only if the comoving number 

density is constant in position and where /i = k-r/k, r be- 
ing the unit vector along the line of sight and k the wave 
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For each redshift bin 










6 


shot noise 


Ps 
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angular diameter distance 


log Da 
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Hubble parameter 


logH 
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growth factor 


logD 


10 


bias 


log/3 


Table I: Cosmological parameters 




Parameters 
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total matter density 


UJm — ^rriQ^ 


2 


total baryon density 


LJb ~ ilbh^ 


3 


optical thickeness 


T 


4 


spectral index 


Us 


5 


matter density today 
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tensor scalar ratio 


T/S 


7 


angular diameter distance 


log Da 


8 


normalization factor 


log As 



Table II: CMB parameters 



vector. The highest frequency kmax{z) is chosen to be 
near the scale of non-linearity at z: we choose values from 
O.ll/i/Mpc for small z bins to 0.33/i/Mpc for the highest 
redshift bins. Any submatrix of F~j^ gives the correla- 
tion matrix for the parameters corresponding to rows and 
columns on that submatrix. The eigenvectors and eigen- 
values of this correlation matrix give the orientation and 
the size of the semiaxes of the confidence region ellipsoid. 
This automatically marginalizes over the remaining pa- 
rameters. The parameters that we use for evaluating the 
Fisher matrix are shown in Tab. Our fiducial model 
corresponds to the ACDM WMAPSy best-fit parameters 
(Spergel et al. 2006): 17,„o = 0.28, h = 0.73, ^de = 0.72, 
= 0, 17fc/i2 = 0.0223, T = 0.092, = 0.96 and 
T/S ~ and as anticipated 7 — 0.545. Beside the 
BAO from large scale structure, we also employ the CMB 
Fisher matrix, following the method in Eisenstein, Hu & 
Tegmark (1999) and assuming a Planck-like experiment. 
The cosmological parameters we use for CMB are listed 
in Tab. The total Fisher matrix is given simply by 
the addition of the two matrices. 

The derivatives of the spectrum with respect to the cos- 
mological parameters Pi (i.e. a;„j = = ^bo^^ 



T, n.. 



r^nio plus Pg, (3,G, D, H for each redshift bin) are 




Figure 1: Growth index vs the scale factor for a DE model 
with a varying equation of state, w{z) = wo + 'Wiz with wo = 
— 1.5 and wi = 1. The black solid line refers to the solution 
obtained by the differential equation for 7 (Percival 2005). 
The red dashed line refers to the growth index given by eq. 
([7]) and the green dotted line is the growth index with a 7 
factor given by eq. (HZ}- The matter density is given by eq. 



Since we want to propagate the errors to the cosmo- 
logically relevant set of parameters 



(13) 



we need to change parameter space. This will be done 
taking the inverse of the Fisher Matrix F^^ and then 
extracting a submatrix, called containing only the 
rows and columns with the parameters that depend on 
qi, namely Da, H and G. The root mean square of the 
diagonal elements of the inverse of the submatrix give the 
errors on Da, H, and G. Then we contract the inverse 
of the submatrix with the new set of parameters qi ; the 
new Fisher matrix will be given by 



F 



DE-ij 



~a — ^ 

oqi 



dpn 

dqj 



(14) 



This automatically marginalizes over all the remaining 
parameters. 

The derivatives of the Hubble parameter and for the 
angular diameter distance can be written as 



i91og H 
dqi 



1 dH 



d\og Da 
dqi 



{1 + z)Da 



d\og H 1 



dz 



(15) 



(16) 



evaluated using the fit of Eisenstein & Hu (1999). 



IV. GROWTH FACTOR 

We consider now separately two cases: in Case 1 the 
growth rate depends on w (assumed constant); in Case 
2 the growth rate is free and we forecast the constraints 
that future experiments can put on it. 
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Surveys 


z 


Vs (Gpc/h)3 


n 


0-0.5 


0.006 


5 • 10"^ 


0.5-0.7 


0.0082 


6.9 • 10-2 


0.7-0.9 


0.011 


4.2 • 10-2 


0.9 - 1.1 


0.0135 


3.1 • 10-2 


1.1 - 1.3 


0.015 


2.4 • 10-2 


2.7-3.5 


0.073 


2 ■ 10-^ 



Table III: Details of the surveys. 



z = l+3; A=1000deg- 




z = l+3; A=5000deg^ 




z = 3; A=1000 deg- 




z = 3; A=5000 deg' 



-1.06 -1.03 -1 -0.97 -0.94 




Figure 3: Fom for Wp — w\ vs marginalized parameters. 





1000 deg^ 


5000 dec/2 


10000 dec/2 




Wp 


Wi 


Wp 


Wi 


Wp 


Wi 




2 = 1+3 


0% 


0.036 


0.313 


0.018 


0.208 


0.012 


0.175 


4% 


0.054 


1.523 


0.028 


0.726 


0.018 


0.51 








z - 


= 3 






0% 


0.044 


0.579 


0.033 


0.257 


0.030 


0.154 


4% 


0.061 


2.614 


0.051 


1.463 


0.034 


1.092 



Table IV: Values of for spectroscopic surveys 

Sz = and photometric surveys 5z = 4% on the redshift 
estimate and for several survey areas (case 1). We consider 
two different combinations of redshift {z = 1 + 3 and z = 3 
only). 



-1.05 -1.025 -1 -0.975 -0.95 



z = l+3;A=10000deg' 




z = 3; A=10000deg^ 




-1.05 -1.025 -1 -0.975 -0.95 



Figure 2: Confidence level for Wp and wi for surveys of 1000, 
5000 and 10000 deg^ and for different combinations of red- 
shift bins (case 1). The solid red curve refers to spectroscopic 
surveys and the dashed blue curve to photometric surveys, 
Sz = and 5z/z = 0.04 respectively. 



Case 1 



In general, the exponent 7 depends on the cosmological 
parameters. To see this, we just need to consider the 
equation of perturbations and insert the growth index 
defined by eq. HI). Then we obtain the approximate 
analytic solution (Wang & Steinhardt 1998) : 



7(z) 



l-w(z) 



(17) 



and for a ACDM model 7 — 0.545. The behavior of the 
growth index for a w{z) model is shown in Fig. |(l|). We 
can see that there is almost no difference in behavior be- 
tween the curves obtained with the approximation (fT7|) . 
Because of the dependence of 7 on the dark energy pa- 
rameters, the derivatives of the growth factor are given 
by: 



d\og G 
dqi 



— log n^{z) 



+7 



dqi 



dz 



(1 



(18) 



In this case the new set of parameters is qi = {wq, wi} 
and we assume as fiducial model wq = —l,wi = 0. The 
factor 7, in this case, depends only on the dark energy 
parameters wq and wi ; this means the only non- vanishing 
derivatives are and S-^- In Fig. |l2| the confidence 

OWQ awi o H — y 

regions are shown for different combination of redshift 
and area. Instead of (wo, wi) we use the pivot parameters 
Wp — wi (projection of Wq- on the pivot point, defined 
as the value of z for which the uncertainty in w (z) is 
smallest). 
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z = 1+3; A = lOOOdeg^ z = 3; A = lOOOdeg^ 




-1.05 -1 -0.95 -1-05 -1 -0.95 

Wo Wq 



z = 1+3; A = 5000deg- z = 3; A = 5000deg^ 




w„ Wo 



z = 1+3; A = lOOOOdeg- z = 3; A = lOOOOdeg' 






1000 deg'^ 


5000 deg^ 


10000 deg'^ 


5z 


Wo 


7 


wo 


7 


wo 


7 




2 = 1+3 


0% 


0.045 


0.099 


0.016 


0.059 


0.004 


0.05 


4% 


0.128 


0.301 


0.062 


0.153 


0.044 


0.114 








z - 


= 3 






0% 


0.089 


0.188 


0.039 


0.092 


0.026 


0.069 


4% 


0.152 


0.344 


0.076 


0.18 


0.081 


0.197 



Table V: Values of o'lUQ and (T^ for spectroscopic surveys 5z — 
and photometric surveys Sz/z = 4% on the measure of 
the redshift and for several areas (case 2). We consider two 
different combinations of redshift bins (z = 1 + 3 and z = 3 
only). 



z=0+l.S; 20000 deg" 



-0.05 
-0.1 




-1.06 -1.04 -1.02 -1 -0.98 -0.96 -0.94 



Figure 4: Confidence level for wq and 7 for surveys of 1000, 
5000 and 10000 deg^ and for different combinations of red- 
shift bins (case 2). The solid red curve refers to spectroscopic 
surveys and the dashed blue curve to photometric surveys, 
Sz = and Sz/z = 0.04 respectively. 

B. Case 2 

We want now to put constraints on 7 as a free param- 
eter. We assume here w — constant and again wq — —1 
as fiducial value. The new set of parameters is there- 
fore Qi = {wo, 7}. The derivatives with respect to the 
first three parameters are given by the eq. ifTSj) . The 
derivative for the growth factor with respect to 7 is: 

d\ogG f d r 1 o M 

= - y log (z) {zf (19) 



V. RESULTS AND CONCLUSIONS 

The main aim of this work is to give marginalized con- 
straints on the dark energy parameters {wp — wi) and 



Figure 5: Confidence level for Wp and w\ for surveys 20000 
deg'^ (DETF case). The solid red curve refers to spectroscopic 
surveys. 



most importantly on the growth factor itself, for several 
combinations of surveys, redshift errors and area. Fol- 
lowing SE and Amendola, Quercellini, Giallongo (2004) 
we consider several binned surveys with average redshift 
depth around z = 1 and z — 3 plus a SDSS-like survey 
at z < 0.5, as detailed in Table Hill More details can be 
found in AQG. We consider both spectroscopic surveys 
{Sz = ) and photometric surveys {Sz/z = 0.04) and 
three areas (1000, 5000, 10000 deg^). These features are 
well within the range of proposed experiments like JDEM 
and DUNE (Crotts at al. 2005; Refregier et al. 2006; see 
also DETF Report Albrecht et al. 2006) 

We first consider Case I, in which the growth factor is 
not an independent quantity but is a function of w{z). 
The two-dimensional regions of confidence are shown in 
Fig |[2| and the final marginalized errors are summarized 
in Tab. IjlVp . The errors on Wp reduce from 0.036 to 
0.012 for the spectroscopic case for surveys that extend 
from 1000 to 10000 deg^ and from 0.054 to 0.018 in the 
photometric case. 

Then we consider Case II, in which 7 is a free constant 
as in eq. ([1]). In Fig. ^ we show the confidence regions 



z=0+1.5; 20000 deg^ 

0.65 , : — 




OAS 

-1.02 -1 -0.98 -0.96 

Wo 

Figure 6: Confidence level for wo and 7 for surveys 20000 
deg'^ (DETF case). The solid red curve refers to spectroscopic 
surveys. 



for Wo — 7- The errors are given in Tab. ([Vl . We see 
that the errors on 7 reduce from 0.099 to 0.05 for spec- 
troscopic surveys and from 0.301 to 0.114 in the photo-z 
case. These errors are way too large to produce an inde- 
pendent constraint on w (in fact one has approximately 
Aw; « IOA7 near w = —1) but, beside being a gen- 
eral test of consistency for the cosmological model, they 
would certainly give interesting constraints on models 
that predict growths different from standard, like modi- 
fied gravity models (Koyama & Maartens 2005; Maartens 
2006; Amendola, Charmousis & Davis 2005; Amendola, 
Polarski, Tsujikawa 2006). In Fig. ^ we show the 
FOM for Wp — wi, for only one survey (5000 deg^) and 
only one combination of redshift (z = 1 — 3), first when 
all the other parameters are fixed and then successively 
marginalizing over the parameter indicated and over all 
those on the left (eg the third column represents the 
marginalization over ujm, ujb). 

We can compare our results to those obtained recently 
by Huterer and Linder (2006). Using a combination 
of weak lensing, SNIa and CMB methods, they predict 
a (7) = 0.044 for future experiments. With large-scale 
tomographic weak lensing alone, Amendola, Kunz, and 
Sapone (2007) predict a (7) = 0.04 at 68% confidence 
level. These values are comparable to those obtained 
here with the BAO method and considering a spectro- 
scopic survey of 5000c?eg^, a (7) = 0.059. 

We notice that the difference on the growth index 7 be- 
tween General Relativity and an extradimensional grav- 
ity model (as DGP, where 7 = 0.68, see Linder & Cahn 
2007) is A7 = 0.135; if we compare our results shown 
in Tab. ||V]) we see that the errors on 7 for a photo- 
metric survey are within this range, meaning that DGP 
model cannot be escluded. Things get slightly better if 
we consider spectroscopic surveys, where errors decrease 
with about 30%; however in this case we require a large 
survey extended from z = 1 to z = 3 with an area of 
10000 de^^ to distinguish with sufficient confidence DGP 
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Figure 7: Effect of Survey geometry on the dark energy FOM. 
We plotter the FOM (wp—wi) for spectroscopic surveys {Sz = 
0) and photometric surveys {Sz = 0.04) as a function of the 
area. 
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Figure 8: Effect of Survey geometry on the dark energy FOM. 
We plotter the FOM {wo —7) for spectroscopic surveys {Sz — 
0) and photometric surveys {Sz = 0.04) as a function of the 
area. 



from ACDM. In Fig. ^ is shown the confidence region 
for Wo — 7 for a survey extended from z = 0toz=1.5 
and an area of 20000 deg^ (DETF case) : the error on 7 
reduces to a (7) = 0.06. 

In Fig. ^ we also show the figure-of-merit (FOM) sug- 
gested by the Dark Energy Task Force report ( Albrecht 
et al. 2006) as a simple measure of the constraining power 
of an experiment. The FOM is defined as the inverse of 
the area that encloses the 95 % confi dence region and can 
be found simply as (6.177r-\/det F)~^. In Fig. ^ we 
plot the FOM for wq and 7. The general trend is that 
the FOM for spectroscopic surveys are roughly 4-6 times 
higher than for similar 4% error photo-z surveys. It will 
be interesting to compare our FOM on the plane ^0,7 
with those obtained from other experiments. This task 
will be performed in future work. 
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